两组率的Meta分析，手把手教Stata操作！

一、问题与数据

Franik S等人拟探讨芳香化酶抑制剂治疗多囊卵巢综合征患者，其临床活产儿率是否高于使用克罗米酚柠檬酸盐治疗的患者。研究者通过检索数据库，对最终纳入的文献，提取了如下信息：作者、发表年代、试验组活产儿数量、试验组样本量、对照组活产儿数量、对照组样本量。

最终纳入6篇研究，合计1353名妇女。各研究提取得到的数据如表1。试问：对于多囊卵巢综合征，使用芳香化酶抑制剂治疗的患者，其临床活产儿率是否高于使用克罗米酚柠檬酸盐治疗的患者？

表1 文献信息提取表

二、对数据结构的分析

该Meta分析纳入了6个研究，每个研究均给出了试验组活产儿数量、试验组样本量、对照组活产儿数量和对照组样本量4个数据。Meta分析的结局指标为两组活产儿率比较的OR值，因此属于两组率的Meta分析。

三、Stata操作方法

1. 创建数据集并录入数据

(1) 打开数据录入窗口：Data → Data Editor → Data Editor (Edit)

关于数据录入细节，可参考文章：Stata做Meta分析：以诊断试验准确性为例 

(2) 录入数据后，点击右上角关闭按钮

各变量名代表含义如下表所示：

（注：非活产儿数量=样本量-活产儿数量）

2. 两组率的Meta分析

(1) 操作方法：在Command窗口中进行编程。

在“Command”栏中输入编程语句→点击“Enter”键

注：在进行分析前，需要先在stata中安装meta分析模块。

在Command栏中输入语句：

ssc install metan

(2) 基本原理

进行两组率的meta分析，实际上为对原始研究中试验组和对照组活产儿率的比值进行合并。

(3) 具体编程步骤：使用固定效应模型进行两组率的meta分析。

程序为：metan aevent anonevent bevent bnonevent, or fixed

程序解读：

两组率的meta分析的程序为：

Metan 试验组发生数 试验组未发生数 对照组发生数 对照组未发生数, 效应指标 使用的模型

metan代表对数据进行meta分析；or代表选择的效应指标为比值比（病例对照比较研究中，两组率的meta分析的效应指标均为or）；fixed代表固定效应模型；（若要选择随机效应模型，应使用random）。

四、结果解读

1. 首先看异质性检验（Heterogeneity）的结果：I²=0%，说明6个原始研究的数据之间不存在异致性，应选用固定效应模型。

2. 再看合并的OR值其95%CI（Model Results）：合并的OR=1.805，95%CI：1.396 -2.333。

下图所示为Meta分析的森林图。结果解读同上。

五、撰写结论

固定效应模型Meta分析结果显示：对于多囊卵巢综合征，使用芳香化酶抑制剂治疗的患者，其临床活产儿率高于使用克罗米酚柠檬酸盐治疗的患者，OR值为1. 805（95%CI：1.396 -2.333）。

六、延伸阅读

1. 模型选择的标准

若各原始研究间存在异质性，则使用随机效应模型；否则，则使用固定效应模型。

2. 异质性的判断

Meta分析时，若I²<25%，则认为不存在异致性；若I²介于25%-50%之间，则认为异致性程度较小；若I²的值介于50%-75%之间，则认为存在一定的异致性；若I²>75%，则认为存在较大的异致性。

本文案例文献来源：

Cochrane Database Syst Rev. 2014, Issule 2. CD010287. 

相关阅读

医咖会微信：medieco-ykh

关注医咖会，轻松学习研究方法学~

有临床研究设计或统计学方面的难题？快加小咖个人微信（xys2016ykf），拉你进统计讨论群和其他小伙伴们一起交流学习；或者点击公众号下方自定义菜单的“统计咨询”，提出你遇到的统计难题。

点击左下角“阅读原文”，看看既往小伙伴们都提出了哪些问题，以及该如何去解决，也许正好能帮到你~